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Abstract 

Recent results on solar neutrinos provide hints that the LMA MSW solution 
could be correct. We perform accurate calculations for potential 'smoking- 
gun' effects of the LMA solution in the SuperKamiokande solar neutrino ex- 
periment, including: (1) an almost constant reduction of the standard recoil 
electron energy spectrum (with a weak, < 10%, relative increase below 6.5 
MeV); (2) an integrated difference in day- night rates (2% to 14%); (3) an ap- 
proximately constant zenith-angle dependence of the nighttime event rate; (4) 
a new test for the difference in the shape of the equally-normalized day-night 
energy spectra (~ 1%); and (5) annual variations of the signal due to the 
regeneration effect (~ 6 times smaller than the integrated day-night effect). 
We also establish a relation between the integrated day-night asymmetry and 
the seasonal asymmetry due to LMA regeneration. As a cautionary example, 
we simulate the effect of an absolute energy calibration error on the shape 
(distortion) of the recoil energy spectrum. 

We compare LMA predictions with the available SuperKamiokande data 
and discuss the possibilities for distinguishing experimentally between LMA 
and vacuum oscillations. If LMA is correct, global solutions combining data 
from different types of measurements made by SuperKamiokande or by dif- 
ferent solar neutrino experiments could reveal in the next few years a many 
a indication of neutrino oscillations. 







I. INTRODUCTION 



On April 1, 1996, solar neutrino research shifted from the pioneering phase of the chlorine 
(Homestake) experiment [l| and the exploratory phase of the Kamiokande 0, SAGE ||, 
and GALLEX [|J experiments, to the era of precision measurements inaugurated by the 
SuperKamiokande [|5HTT|1 experiment. In this paper, we concentrate on the predictions for 
'smoking-gun' evidence of new neutrino physics that can be observed in electron-neutrino 
scattering experiments like SuperKamiokande. 

There are four generic oscillation solutions (LMA, SMA, and LOW MSW solutions and 
vacuum oscillations, cf. Ref. |T^] for a recent discussion) of the solar neutrino problems 
that involve two neutrinos. These solutions are all consistent with the predictions of the 
standard solar model |]T3[ and the observed average event rates in the chlorine, Kamiokande, 
SuperKamiokande, GALLEX, and SAGE experiments. We limit ourselves to such globally 
consistent oscillation solutions and adopt the methodology used in our recent paper |T2| . 



Recent experimental results from SuperKamiokande |6|-|9] provide some encouragement 
for considering the LMA solution of the MSW effect ||14|| . No statistically significant dis- 
tortion of the recoil electron energy spectrum has been discovered. No obvious change 
has been revealed in the slope of the ratio [Measurement/ (Standard Spectrum)] as the 
SuperKamiokande measurements have been extended to progressively lower energies. More- 
over, the excess of events at the highest energies has diminished. 

A larger flux has been observed at night than during the day, although the difference 
observed by SuperKamiokande @-|| (which is about the size of what is expected from the 
LMA solution) is only ~ 1.6a after 708 days of data collection ||15|| . In addition, no large 



enhancement of the rate has been found when the neutrino trajectories cross the earth's 
core. Moreover, the excess that is found appears already when the sun is just below the 
horizon. These results are, as we shall see in detail in this paper, expected on the basis of 
the LMA MSW solution but are not expected on the basis of the SMA solution. 

There are also some non-solar indications that large mixing angles may be plausible. 
First, and most important, SuperKamiokande and other experiments on atmospheric neu- 
trinos have provided strong evidence that large lepton mixing is occurring at least in one 
channel involving muon neutrinos fL6 |. Moreover, the mass required for the LMA solution 



(m ~ y (Am 2 ) ~ 10~ 2 eV) is only one order of magnitude smaller than the mass scale sug- 
gested by atmospheric neutrino oscillations. All of the other solar neutrino solutions imply 
a smaller or much smaller neutrino mass scale for solar neutrinos. It is however well known 
that a weak mass hierarchy allows a rather natural explanation of large mixing. 

The oscillation parameters of the LMA solar neutrino solution can give an observable ef- 



fect for atmospheric neutrino experiments. In particular, the LMA solution could explain (T? 
the excess of e— like events in the atmospheric data. 

In this paper, we provide accurate calculations for smoking-gun predictions of the LMA 
solution, including a simple new test (Sec. [Ill Q ). In particular, we calculate and compare 
with SuperKamiokande observations the following quantities: (1) the distortion of the elec- 
tron recoil energy spectrum; (2) the zenith angle dependence of the observed counting rate; 
(3) the day-night spectrum test; and (4) the seasonal dependence (beyond the effect of the 
earth's eccentric orbit) of the counting rate. All of these quantities are zero, independent of 
solar physics, if the standard model for electroweak interactions is correct. Therefore, mea- 
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TABLE I. Standard Solar Model and LMA Predictions for Solar Neutrino Rates. The standard 
solar model predictions are taken from BP98, Ref. |f[3[j . The LMA predictions refer to the best-fit 
LMA solution when only the total rates are considered: Am 2 = 2 x 10 _5 eV 2 and sin 2 26 = 0.8 
(cf. Ref. p^l). [The values given in parentheses were computed using, respectively, the oscillation 



parameters: Am 2 = 8 x 1(T 5 eV 2 ,sin 2 26» = 0.63 and Am 2 = 8 x 10" 5 eV 2 ,sin 2 (9 = 0.91.] The 
rates calculated using the BP98 fluxes and the LMA oscillation solution agree with the measured 
rates for the Homestake chlorine experiment (2.56 ± 0.23 SNU, see Ref. Q) and the SAGE (67 ± 8 
SNU, see Ref. [§) and GALLEX (78 ± 6 SNU, see Ref. §) gallium experiments. 



Source 


Flux a 
Standard 


CI 

(SNU) 
Standard 


Ga 
(SNU) 
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CI 
(SNU) 
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Ga 
(SNU) 
LMA 


PP 

pep 
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7 Be 

8 B 

13 N 
15Q 


5.94 

1.39 x 10~ 2 
2.10 x 10~ 7 
4.80 x 10" 1 
5.15 x 10~ 4 
6.05 x 10~ 2 
5.32 x 10~ 2 


0.0 

0.2 

0.0 

1.15 

5.9 

0.1 

0.4 


69.6 
2.8 
0.0 
34.4 
12.4 
3.7 
6.0 


0.0 (0.0,0.0) 
0.1 (0.1,0.1) 
0.0 (0.0,0.0) 
0.46 (0.65,0.60) 

1.8 (1.6,2.3) 
0.05 (0.06,0.05) 
0.16 (0.25,0.20) 


39.8 (47.1,37.5) 
1.2 (1.8,1.4) 
0.0 (0.0,0.0) 

15.7 (21.5,17.4) 

3.7 (3.4,4.8) 

1.8 (2.4,1.9) 
2.5 (3.8,3.1) 


Total 




7 7+1.2 
'•'-1.0 


129+J 


2.6 (2.65,3.2) 


64.7 (80,66) 


a All fluxes 


are in units of 


(10 10 cm-V 1 ) 









surements of the spectral energy distortion and the zenith angle and seasonal dependences 
constitute sensitive smoking-gun tests of new physics. As a cautionary example, we also 
simulate the effect of an energy calibration error on the apparent distortion of the electron 
recoil energy spectrum. 

Figure [I] shows the allowed region of LMA solution in the plane of sin 2 29 and Am 2 . We 
have indicated the region allowed by the measured total rates in solar neutrino experiments 
by the continuous contours in Fig. [IJ the contours are shown at 90%, 95%, and 99% C.L. 

The constraints provided by the total rates in the Homestake, GALLEX, SAGE, and 
SuperKamiokande experiments are, at present, the most robust limitations on the allowed 
parameter space. The only slight change in the available data on total rates since the publica- 



tion of Ref. |L2| is a small reduction after 708 days in the total error for the SuperKamiokande 
experiment ||; the allowed region based upon the total rates in solar neutrino experiments 
is essentially unchanged from our earlier study (cf. the similar results in Fig. 2 of Ref. [12| 
for the case that only the total rates are considered). 

The dashed line in Fig. |I| represents the 99% allowed contour including both the measured 
SuperKamiokande day-night effect and the total rates. The SuperKamiokande collaboration 
has performed preliminary global analyses of the available 708 days of data PJlO,15 . 



Table [I] contrasts the predictions of the standard model and the LMA for the chlorine 
and the gallium experiments. The LMA results that are not in parentheses are for Am 2 = 
2 x 10~ 5 eV 2 and sin 2 29 = 0.8. These neutrino parameters are close to the best-fit LMA values 
when only the total rates of the chlorine, Kamiokande, SuperKamiokande, GALLEX, and 
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SAGE experiments are considered (see Ref. [fL2] ). In parentheses, we show the predictions 
for two rather extreme LMA solutions. We have evaluated the predictions for the chlorine 
and gallium experiments for a representative set of oscillation parameters and all of the 
values lie within the boundaries defined by the three solutions shown in Table |. Therefore, 
the predictions given in the table show the expected range of capture rates consistent with 
the LMA solution and the existing experiments. The range of production rates predicted by 
the set of LMA oscillation solutions considered in this paper is: 37 C1 = 2.9 ± 0.3 SNU and 
71 Ga = 71 ± 9 SNU. 

The prediction of the LMA solution shown in Table | for the SuperKamiokande rate due 
to the 8 B neutrino flux [|,|||| is : Rate (measured, 8 B) = 0.474 Rate (BP98, 8 B). 

This paper is organized as follows. We discuss the distortion of the recoil electron energy 
spectrum in Sec. Til the day-night differences of the event rate in Sec. [Til], and the seasonal 
dependences in Sec. [TV|. We describe in Sec. [V| possibilities for distinguishing between LMA 
and vacuum oscillations. We summarize and discuss our main results in Sec. [VI]. 



II. SPECTRUM DISTORTION 

Figures |2] and ^| illustrate the main results of this section. The reader who wants to get 
quickly the essence of the physical situation can skip directly to the discussion of Fig. ^ in 
Sec. ffg and Fig. | in Sec. TO 

The observation of a distortion in the recoil electron energy spectrum produced by so- 
lar neutrinos scattering off electrons would be a definitive signature of physics beyond the 
Standard Electroweak Model. Neutrinos which correspond to recoil electrons with energies 
between 5 MeV and 13 MeV, which are most easily observed by the SuperKamiokande solar 
neutrino experiment ||[6|||[|, are predicted by the Standard Solar Model to be essentially 
all (more than 99%) from 8 B beta decay [|13j. The shape of the 8 B neutrino energy spectrum 



can be determined accurately from laboratory measurements [31| and the influence of astro- 



physical factors is less than 1 part in 10 5 For energies relevant to solar neutrino studies, 



the cross-sections for neutrino scattering have been calculated accurately in the standard 
electroweak model, including radiative corrections [^(J. Therefore, if nothing happens to 
solar neutrinos on the way from the region of production in the solar interior to a detector 
on earth, the recoil electron energy spectrum from 8 B neutrinos can be calculated precisely. 

Any measured deviation of the recoil electron energy spectrum from the standard shape 
would indicate new physics. The opposite, however, is not true: the absence of a measured 
distortion does not mean the absence of new physics. 

We first indicate in Sec. [11 A| how the neutrino energy spectrum is convolved with the 



electron recoil energy spectrum and with the measurement characteristics of the detector. 



We then summarize briefly in Sec. [IIB| the characteristics of the measured SuperKamiokande 
energy spectrum. In Sec. [11 Q we answer the question: How does an error in the energy 
calibration affect the shape of the recoil energy spectrum? We discuss in Sec. [II D| the 
comparison of the predicted and the measured electron energy spectra, which are presented 
in Fig. H. 
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A. The convolution 



We concentrate in this paper on analyzing the ratio, R(Ei), of the number of observed 
events, Np hs in a given energy bin, E^ to the number, Nf SM , expected from the SSM [13] , 
where 

N° hs 

For specific applications, we use the observed event numbers in the form provided by Su- 
per Kamiokande ||||; there are 17 energy bins of 0.5 MeV width from 5.5 MeV to 14 MeV 
and an 18th bin that includes all events from 14 MeV to 20 MeV. When more data are 
available, it will be important to divide the events above 14 MeV into several different 



bins [21 



The number of events in a given energy bin, i, can be expressed as 

rEi+O.BMeV 



rEi+0.5MeV roo roo 

K= dE e dE'J(E e ,E' e ) dE u F(E u )x 

Je z Jo J el 



^ da e {E u ,E' e ) dap(E v ,E' e ) 
P[Eu} dE< + (1 " P[Eu)) dE' 



(2) 



where F(E V ) is the flux of 8 B neutrinos per unit energy at the detector, f(E e ,E' e ) is the 
energy resolution function which we take from KJ. P{E U ) is the survival probability v e — > 
v e , da e /dE' e and da^/dE' e are the differential cross-sections of the v e e— and v^e— scattering. 
There are practically no 8 B neutrinos with energies greater than 15.5 MeV, so that for 
practical purposes the upper limit of the integral over E v in Eq. (0) can be taken to be 16 
MeV. Including the effects of the finite energy resolution, the upper limit over the measured 
electron recoil energy, E' e , can be taken to be 20 MeV. 

To use Eq. (0) to calculate the combined prediction for the standard electroweak model 
and the standard solar model, set P = 1. The second term in the brackets of Eq. (0) 
then disappears. This term is also absent if the electron neutrinos are converted to sterile 
neutrinos. 

How does the convolution shown in Eq. (Q) affect energy-dependent distortions (or con- 
versions) of the incoming electron neutrinos? The convolution spreads out over a relatively 
wide energy range any energy- dependent features. In particular, the scattering process pro- 
duces electrons with energies from zero to almost equal to the incoming neutrino energy (for 
neutrino energies much greater than the electron mass). If v e are converted to or /and u T , 
then there will be neutral current scattering [the second term in boxed brackets of Eq. @] 
which, while less probable, will also reduce the apparent effect of the conversion. Since the 
solar 8 B neutrino spectrum decreases rapidly beyond about 10 MeV, a distortion of the recoil 
electron spectrum at some energy E e (E e > 10 MeV) is determined by the distortion of the 
neutrino spectrum at somewhat lower energy E v < E e . The energy resolution function has 
a crucial affect in smearing out distortion. The 2a width for the energy resolution is about 
3 MeV at E e = 10 MeV in SuperKamiokande [§,H§. Smaller features will be smeared out. 
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B. Summary of the data 



After 708 days of data taking with SuperKamiokande, no unequivocal distortion of the 
electron recoil energy spectrum has been found. The data show essentially the spectrum 
shape expected for no distortion for E < 13 MeV with some excess events at E > 13 MeV. 
The excess electrons observed at higher energies may reflect: (1) a statistical fluctuation; (2) 
the contribution of hep neutrinos [PH -231; (3) a larger-than-expected error in the absolute 
energy normalization; or (4) new physics. The systematic effects due to hep neutrinos or 
to an error in the energy normalization can hide a distortion due to neutrino conversion 
or cause an artificial distortion. For example, a suppression due to oscillations of the high 
energy part of the electron spectrum could be compensated for (hidden by) the effect of hep 
neutrinos. 

All of the four explanations for the high-energy excess have been described previously in 
the literature. However, the distortion due to an error in the absolute energy scale of the 
electrons has only been mentioned as a possibility |[24|| . Therefore, we discuss in Sec. |11 C 



the effect on the energy spectrum of an error in the absolute energy scale. 

Another possible systematic error, a non-Gaussian tail to the energy resolution, could 
also produce an apparent distortion at the highest energies. However, this effect has been 
estimated by the SuperKamiokande collaboration and found to be small 

To resolve the ambiguity at high recoil electron energies, one can proceed in at least two 
different ways. (1) One can exclude the spectral data above 13 MeV, since the distortion due 
to both hep neutrinos and the calibration of the energy scale normalization occur mainly at 
high energies. A comparison of the results of the analysis using all the data with the infer- 
ences reached using only data below 13 MeV provides an estimate of the possible influence 
of systematic uncertainties. However, this procedure does throw away some important data. 
(2) One can treat the solar hep neutrino flux as a free parameter in analyzing the energy 



spectrum |2T|j22]| . In general, this is the preferred method of analysis and the one which we 
follow. 



C. How does an energy calibration error affect the recoil energy spectrum? 

An error in the absolute normalization of the energy of the recoil electrons, 5 = E trne — 
-^measured, would lead to an apparent energy dependence in R that is particularly important 
for the LMA solution, since the theoretical prediction is that R is essentially constant in the 
higher energy region accessible to SuperKamiokande. The ratio R{E) can be written as 

R obs m N°*(E + 6) „ t / 1 dN SM (E) \ 

R {E) ~ N™(E) ~ R ° {E) y + N™(E) d E 5 ) ' (3) 

where N SM (E) is the standard model spectrum. Since N SM , dN SM ^ dE and, in general, 
6 all depend on energy, the last term in Eq. (|3|) gives rise to an apparent distortion of 
the observable spectrum. Moreover, since N SM (E) is a decreasing function of energy, 5 x 
(dN SM /dE)/N SM increases with energy for negative 5 . Thus, for negative 5, R{E) is 
enhanced at high energies. The SuperKamiokande collaboration has estimated the error in 
the absolute energy calibration as being 0.8% at 10 MeV M, i.e., 5 = 80 keV, la. 
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We have simulated the effect of an error in the absolute energy calibration by con- 
volving a standard model recoil energy spectrum with the energy resolution function of 
Super Kamiokande and introducing a constant offset error S. We calculated the ratio R(E) 
directly from the equality in Eq. (0), without making the approximation involved in the 
Taylor-series expansion. 

Figure || compares the computed R(E) for three different values of 5 with the observed Su- 
perKamiokande data M for 708 days. For visual convenience, each curve is normalized some- 
what differently, but all are normalized within the range allowed by the SuperKamiokande 
measurement of the total rate. We see from Fig. [2] that a 2a (160 keV) energy calibration 
error can produce a significant pseudo-slope. More exotic effects can be produced if the one 
assumes that the energy error, 5, is itself a function of energy. 

We conclude that future analyses of the distortion of the energy spectrum should include 
5 as one of the parameters that is allowed to vary in determining from the measurements 
the best-fit and the uncertainty in the energy distortion. 

D. Predicted versus measured spectra 

Figure ^| compares the predicted versus the measured ratio R(E) [see Eq. ([!])] of the 
electron recoil energy spectra for a representative range of Am 2 that spans the domain 
permitted by the global LMA solutions. All of the solutions shown have sin 2 26 = 0.8 and 
are consistent with the average event rates in the chlorine, Kamiokande, SAGE, GALLEX, 
and SuperKamiokande experiments. Also, all of the histograms of neutrino solutions are 
normalized so as to yield the same value of the ratio R(E) at E = 10 MeV. The measured 
values are taken from ||. 

The most conspicuous feature of Fig. |3] is the flatness of the ratio R(E) that is predicted 
by the LMA solution with no enhancement of the hep flux. The survival probability, P, 
without earth regeneration is P(E) ~ sin 2 6 in the observable high energy part of the 8 B 
neutrino spectrum. Regeneration in the earth can lead to a weak increase of the survival 
probability with energy. 

Excellent fits to the recoil energy spectra can be obtained by increasing the production 
cross section for the hep flux by a factor of 10 to 40 over its nominal (see Ref. |p6| ) best 
value. This result is illustrated in Fig. ||. 

Results similar to Fig. |3] are obtained if sin 2 26 is allowed to vary over a representative 
range of the allowed global LMA solutions that are consistent with the average measured 
event rates. We find numerically that the dependence of the shape of the predicted spectrum 
upon sin 2 26 is very weak. 

Recent work by the SuperKamiokande collaboration has placed an experimental upper 
limit on the hep neutrino flux at earth of at most 20 times the nominal standard value |||9|,|l0|] . 
The fits shown in Fig. |3| demonstrate that we require a hep flux at the sun of 10 to 40 times 
the standard value. However, LMA neutrino oscillations reduce the rate at the earth so that 
the observed rate is about a factor of two less than would correspond to a flux of purely v e 
hep neutrinos. So the fits shown for variable hep flux in Fig. |3] and are consistent with the 
SuperKamiokande upper limit on the measured hep flux. 

What are the chances that SuperKamiokande can obtain smoking-gun evidence for a 
distortion of the energy spectrum if the LMA solution is correct? The lack of knowledge of 
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the hep production cross section prevents fundamental conclusions based upon the suggestive 
higher energy (> 13 MeV) events. The data reported by SuperKamiokande (at the time 
this paper is being written) above 14 MeV are in a single bin (14 MeV to 20 MeV). As 
emphasized in Ref. |[21|| , it is possible that measurements in smaller energy bins above 14 
MeV could reduce the uncertainty in the hep flux and test the consistency of the LMA plus 
enhanced hep predictions (for preliminary results see Ref. U). 

For larger allowed values of Am 2 , the distortion curve turns up at low energies (cf. 
Fig. |3|). This upturn is due to the fact that for larger Am 2 the low energy part of the 8 B 
neutrino spectrum (E ~ 5 — 6 MeV) is on the adiabatic edge of the suppression pit. For the 
smallest allowed Am 2 , LMA predicts a weak positive slope for R(E) that is caused by the 
Earth regeneration effect. 

The predicted distortion at low energies, between, e.g., 5 MeV and 6.5 MeV, is small 
(typically ~ a few percent up to as large as 10 percent, cf. Fig. ^) and would require 
approximately 10 years of SuperKamiokande data in order to show up in a clear way. 



III. EARTH REGENERATION EFFECTS 



If MSW conversions occur, the sun is predicted to be brighter at night in neutrinos than 
it is during the day ]27| at those energies at which the survival probability of v e is less than 
one-half. This phenomenon, if observed, would be a dramatic smoking-gun indication of 
new physics independent of solar models[|. 

When the sun is below the horizon, neutrinos must traverse part of the earth in order to 
reach the detector. By interacting with terrestrial electrons, the more difficult to detect or 
v T can be reconverted to the more easily detectable v e . The opposite process can also occur. 
The physics of this 'day-night' effect is well understood p7| . We adopt the methodology of 
Ref. ||28|| . The SuperKamiokande collaboration has placed constraints on neutrino oscillation 
parameters from the measurement of the zenith angle dependence in 504 days of data |§ . 

We first discuss in Sec. [Ill A| the day-night asymmetry averaged over all zenith angles 



and integrated over all energies and seasons of the year. We then analyze in Sec. |! IB] the 
dependence of the total rate upon the solar zenith angle. We then describe in Sec. |111 C| a 
new test, which we call the 'day-night' spectrum test. 



A. The Integrated Day-Night Effect 



The average day-night asymmetry measured by the SuperKamiokande collaboration 
is p| 



A„- d 



night — day 
night + day 



0.060 ± 0.036(stat) ± 0.008(syst). 



(4) 



*If neutrino oscillations occur, the predicted zenith-angle dependence changes slightly for different 
solar models since the flavor content of the calculated solar flux arriving at earth is influenced 
somewhat by the assumed standard model fluxes. See Ref. [28|, Sec. XB. 
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Here 'night' ('day') is the nighttime (daytime) signal averaged over energies above 6.5 MeV 
and averaged over all zenith angles and seasons of the year.Q 

Table |I| shows the calculated values for the day-night asymmetry for a range of different 
LMA solutions. For solutions within the allowed LMA domain, we find A n _ d = 0.02 to 0.14 

The dependence of the asymmetry on Am 2 can be approximated with the LMA solution 
space by 



A n - d = 0.2 



10~ 5 eV 2 
Am 2 



(5) 



for sin 2 29 = 0.8 and for the range of Am 2 shown in Table [HJ. This approximation depends 
only weakly on mixing angle. The la interval of A n ^ d that follows from Eq. (|4|) leads to an 
allowed range of 



Am 2 = (2 - 8) x 10" 5 eV 2 , la . (6) 

Figure [I] shows the approximately horizontal lines of equal Night-Day asymmetry in a 
Am 2 - sin 2 29 plot together with region in parameter space (dashed contour) that is allowed 
by a combined fit of the total rates and the Night-Day asymmetry. The best-fit oscillation 
parameters for the combined fit are 

Am 2 = 2.7 x 10~ 5 eV 2 ; sin 2 29 = 0.76. (7) 

The Night-Day asymmetry is 8% for these best fit parameters. 



B. The Zenith Angle Dependence 

Figures f| and || compare the predicted and the observed (after 708 days) dependence of 
the SuperKamiokande || event rate upon the zenith angle of the sun, 9. The results are 
averaged over one year. In order to make the figures look most similar to figures published 
by the SuperKamiokande collaboration, we have constructed the plot using the nadir angle 
6at = 7r — 0. The predicted and the observed rates are averaged over all energies > 6.5 
MeV. 

Three conclusions can be drawn from the results shown in Fig. |4] and Fig. |5|. (1) The 
experimental error bars must be reduced by about a factor of two (total observation time of 
order eight years) before one can make a severe test of the average zenith-angle distribution 
predicted by the LMA solution. (2) Nevertheless, the available experimental results provide 
a hint of an effect: all five of the nighttime rates are larger than the average rate during 



^The definition of A n _ d given in Eq. (Q) is twice as large as in Ref. p8|] ; we changed our definition 
to the one given in the present paper in order to conform to the usage by the SuperKamiokande 
collaboration. Similarly, we do not discuss the valuable constraints provide by moments of the 
day-night effect because the SuperKamiokande collaboration has not yet provided an analysis of 
their results in terms of moments. 
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TABLE II. Day-Night Asymmetries 



Am 2 


A 


At, A 


Ar, A 


(lO" 5 eV 2 ) 


sin 2 2(9 = 0.9 


sin 2 26 = 0.8 


sin 2 2(9 = 0.7 


8 


0.019 


0.018 


0.017 


4 


0.049 


0.050 


0.049 


3 


0.069 


0.073 


0.071 


2 


0.100 


0.107 


0.107 


1.6 


0.126 


0.135 


0.136 



daytime. (3) The zenith-angle dependence predicted by the LMA solution is relatively flat, 
i.e., the flux is approximately the same for all zenith-angles. 

The predicted LMA enhancement begins in the first nighttime bin and the enhancement 
is not significantly increased when the neutrinos pass through the earth's core (cos(G) > 0.8) 
and can even decrease for some parameter choices. As Am 2 decreases, the oscillation length 
increases and the effect of averaging over oscillation phases becomes less effective with the 
result that some structure appears in the calculated zenith angle dependence. For Am 2 < 
lCT 5 eV 2 , the departure from a flat zenith angle dependence becomes relatively large. 

Certain solutions of the SMA, those with sin 2 29 > 0.007, predict a strong enhancement 
in the core that is not observed and therefore these SMA solutions are disfavored B. 

Figures f| and |5] show that for Am 2 = (3 — 4) x 10 -5 eV 2 and a wide range of sin 2 29 the 
LMA solution provides an excellent fit to the zenith angle dependence. 

The theoretical uncertainties |28| (due to the density distribution and composition of the 
earth and the predicted fluxes of the standard solar model) are at the level of 0.2%, about 
an order of magnitude less than the effect that is hinted at by Fig. [| and Fig. |5|. 



C. The Day-Night Spectrum Test 

The spectral dependence of the day-night effect provides additional and independent 
information about neutrino properties that is not contained in the difference between the 
total day and the night event rates. The new physical information that is described by the 
difference between the day and the night energy spectra is the extent to which the different 
energies influence the total rates. 

Figure ^| compares the calculated electron recoil energy spectra for the day and the night 
for a range of values of Am 2 and sin 2 29. Again, we have plotted on the vertical axis of 
Fig. I the ratio, R(E) [see Eq. ([!])], of the measured recoil energy spectrum to the spectrum 
expected if there is no distortion 0. The night spectra lie above the day spectra for all 
relevant energies. Moreover, the difference of the nighttime and the daytime rates increases 



^The average energy spectra during the day and during the night have been calculated also in 
Ref. |2l|. However, these authors did not include the SuperKamiokande energy resolution function, 
which leads to a significant smearing effect (see the discussions in Refs. 
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with energy, reflecting the fact that for the LM A solutions the regeneration effect increases 
with Am 2 . 

The systematic difference between the day and the night spectra is most easily isolated 
when the difference due to the average rates is removed from both the day and the night 
spectra. 

Figure [7| compares the daytime and the nighttime spectra when both spectra are normal- 
ized to the same total rate (the observed SuperKamiokande rate). This equal- normalization 
removes the difference that is normally referred to as the 'day-night effect'. We refer to the 
equal-normalization procedure as the 'day-night spectrum test'. 

Figure [7| shows that for energies below (above) 10 MeV, the predicted differential daytime 
spectrum is higher (lower) than the predicted differential nighttime spectrum. The two 
normalized spectra become equal in rates at about 10 MeV. LMA predicts that the nighttime 
spectrum contains relatively more high energy electrons than the daytime spectrum. 

Figure |7] shows that, when both spectra are normalized to have the same total rates, the 
average daytime spectrum between 5 MeV and 10 MeV is, for the parameters chosen, about 
~ 1.5% greater than the average nighttime spectrum. The oscillation parameters used to 
construct Fig. [7] are Am 2 = 2 x 10~ 5 eV 2 and sin 2 29 = 0.8. We have plotted figures similar 
to Fig. ^ for a number of LMA solutions. For a fixed Am 2 , the equally normalized spectra 
are all very similar to each other. However, the amplitude of the difference between the 
night and the day spectra, about 3% at 5 MeV for the example shown in Fig. [7|, decreases 
to about 2% for Am 2 = 4 x 10" 5 eV 2 and is only about 0.5% for Am 2 = 8 x 10~ 5 eV 2 . 

Does the day-night spectrum test provide information independent of the information 
obtained from the integrated day-night rate effect or the more general zenith-angle effect? 
Yes, one can imagine that the day-night integrated effect (cf. Sec. Ill Aj) has been measured 
to be of order 6% in agreement with LMA predictions, but that when the day-night spectrum 
test is applied the appropriately normalized day rate at energies less than 10 MeV lies below 
the similarly normalized nighttime rate. This later result would be inconsistent with the 
predictions of the LMA solution and therefore would provide information not available by 
just measuring the integrated difference of day and of night rates. 

The easiest way to apply the day-night spectrum test is to normalize the day and night 
spectra to the same total number of events and then compare the number of day events 
below 10 MeV with the number of night events below 10 MeV. In principle, one could divide 
the data into a number of different bins and test for the similarity to the predicted shape 
shown in Fig. |7|, but the Poisson fluctuations would dominate if the data were divided very 
finely. 

The change in slope with energy, illustrated in Fig. [7], between equally normalized day 
and night recoil energy spectra may provide a new test of the LMA solution. It will also be 
useful to calculate the predicted change in slope for other solar neutrino solutions. 



IV. SEASONAL DEPENDENCES 

If the LMA solution is correct, seasonal dependences occur as a result of the same physics 
that gives rise to the zenith-angle dependence of the event rates. At nighttime, oscillations in 
matter can reconvert or v T to v e or matter interactions can also cause the inverse process. 
The seasonal dependence arises primarily because the night is longer in the winter than in the 
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summer. Early discussions of the seasonal effect with applications to radiochemical detectors 
is given in Refs. [31]. A recent discussion of seasonal effects is presented in Ref. [32]]. 



A. Predicted versus observed seasonal dependence 

Figure § shows the predicted LMA variation of the total event rates as a function of 
the season of the year. In constructing Fig. ^ we corrected the counts for the effect of the 
eccentricity of the earth's orbit. We show in Fig. |8| the predictions for a threshold of 6.5 MeV; 
Fig. shows similar results for a threshold of 11.5 MeV. In both cases, the seasonal effects 
are small, ~ 1% to 2%, although the higher-energy events exhibit a somewhat larger (30% 
or 40% larger) variation. The characteristic errors bars for SuperKamiokande measurements 
three years after beginning the operation are larger, typically ~ 4% (cf. Fig. |9]), than the 
predicted LMA seasonal variations. 

How do the SuperKamiokande data compare with the observations? For a threshold 
energy of 6.5 MeV, the data corrected for the earth's eccentricity are not yet available. 
Therefore, we compare in Fig. |S] the observed |9[ and the predicted seasonal plus eccentricity 
dependence. The eccentricity dependence is larger (for the allowed range of parameters) 
than the predicted LMA seasonal dependence; hence, the predicted variations in Fig. |^ 
are larger than in Fig. |[ The regeneration effect enhances the eccentricity effect in the 
northern hemisphere. In the southern hemisphere, the eccentricity and the regeneration- 
seasonal effects have opposite signs. This difference is in principle detectable; if the seasonal 
asymmetry is 7% in the northern hemisphere then it would be about 3% in the southern 
hemisphere. 

Unfortunately, the existing statistical error bars are too large to show evidence of either 
the eccentricity or the LMA seasonal dependence. 

If the LMA solution is the correct description of solar neutrino propagation, it appears 
likely from Fig. |9] that of order 10 years of SuperKamiokande data will be required in order 
to see a highly significant seasonal effect due to LMA neutrino mixing. 



B. The relationship between the summer-winter and the day-night effect 

Since the physics underlying the day-night and the seasonal dependences is the same, i.e., 
non-resonant matter mediated neutrino oscillations, there must be a relationship between 
the two effects. We outline here a brief derivation of a formula that connects the size of the 
day-night asymmetry defined in Eq. and the winter-summer asymmetry, which we define 

as 

_ ^ / Winter — Summer 
U ^ V Winter + Summer 




The principal results given in this section, including Eq. (|To|), have been presented and discussed 



by A. Yu. Smirnov at the Moriond meeting for 1999 |3^]. The seasonal variations found in Ref. [32] 
are in good agreement with our results for the same values of the oscillation parameters. 
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Here 'Winter' ('Summer') is the signal averaged over the period from November 15 to Febru- 
ary 15 (May 15 to August 15). 

Neutrinos reaching the earth from the sun will be in an incoherent mixture of mass 



eigenstates [p4| , p5| . When the sun is below the horizon, the presence of electrons in the 
earth will cause some transitions to occur between different mass states. By solving the 
problem of matter induced neutrino transitions in a constant density medium, one can see 
that the characteristic oscillation length in the earth is always less than or equal to the 
vacuum oscillation length divided by sin 26. The transition probability is proportional to 
sin 2 (</>), where 

Tripsin 20 nn ( R \ /10MeV\ / Am 2 \ 

where R is the distance traversed in the earth. For any of the nighttime bins discussed 



in Sec. |T|, the phase of the oscillation is large so that even relatively small changes in 
distance and energy will lead to fast oscillatory behavior that causes the survival probability 
to average to a constant in all of the bins (cf. Eq. |9|). This is the reason why the zenith-angle 
enhancement shown in Fig. f| is approximately a constant, independent of zenith-angle. 

The event during the night and the event rate during the day are each approximately 
constant and may be represented as Rn and Rd, respectively. Let the average duration of 
night be ts hours during the summer and tw hours during the winter. Then the summer 
signal is proportional to [ii/v^s + -Rd(24 — ts)]- There is a similar expression for the winter 
signal: [R^w + -Rd(24 — tw)]. Then simple algebra shows that the seasonal asymmetry is 

^w-s = (10) 



24 

Thus the seasonal variations are proportional to the night-day asymmetry. The larger the 
day-night effect, the larger is the LMA predicted seasonal variations due to regeneration. 
For the location of SuperKamiokande, tw = 13.9 hr and ts = 10.1 hr, and the parenthetical 
expression in Eq. flTOP is about 1/6. 

Equation ([K]) gives the approximate relation between the seasonal and the day-night 
asymmetries and makes clear why the seasonal dependence is about 6 times smaller than 
the already small day-night effect. For the central value of the measured range of day-night 
asymmetries |1| (cf. Eq. £|), we find from Eq ( |10|) iw-s ~ 1 %, in agreement with our 
explicit calculations. 



V. VACUUM VERSUS LMA OSCILLATIONS 

Table [ITJ summarizes the most striking differences in the predictions of the LMA and the 
vacuum oscillation solutions of the solar neutrino problems. We now discuss some aspects 
of these differences in more detail. 

The sharpest distinctions between vacuum and LMA predictions are expressed in the 
day-night differences. The LMA predicts non-zero day-night differences for all three of the 
tests discussed in Sec. [I ll| and listed as the first three phenomena in Table PTT| . Vacuum 
oscillations predict zero for all these day-night phenomena. 
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TABLE III. LMA versus vacuum oscillations. Section til contains a quantitative discussion of 



the day-night integrated effect, the zenith angle dependence of the rate, and the day-night spectrum 
test. Seasonal effects are discussed in Sec. [V and spectrum distortion is discussed in Sec. ||. 



Phenomenon 



LMA 



Vacuum 



Day-night integrated effect 
Zenith-angle dependence of rate 
Day-night spectrum test 
Spectrum distortion 
Seasonal effect 



small but non-zero 
small but non-zero 
small but non-zero 
flat relative spectrum 
eccentricity dominates 



zero 
zero 
zero 



can explain upturn at large energies 
can be large and energy dependent 



The spectrum distortion is predicted to be small for the LMA solution (cf. Fig. 
The ratio of the observed to the standard spectrum is essentially constant for energies 
above 6.5 MeV, although a modest upturn can occur at lower energies. Vacuum oscillations 
allow a more varied set of spectral distortions and can accommodate, without enhanced hep 
production, the possibly indicated upturn in the relative spectrum beyond 13 MeV. 

One can use seasonal variations to distinguish between vacuum oscillations and LMA 
oscillations. For an early discussion of this possibility, see Ref. The LMA solution 

predicts that the seasonal effects are smaller than the geometrical effect arising from the 
eccentricity of the earth's orbit. Moreover, there is only a weak enhancement of the LMA 
seasonal effect with increasing energy threshold. On the other hand, for vacuum oscillations 
the seasonal effects can be comparable with the geometrical effect and there can be a strong 
enhancement (or an almost complete suppression, even a reversal of the sign) of the seasonal 
effect as the threshold energy is increased (see first paper in Ref. |36| ) . 

For the LMA solution, regeneration occurs but only during the night. Therefore, the 
LMA solution predicts that there is no seasonal dependence of the daytime signal beyond 
that expected from the eccentricity of the earth's orbit. The vacuum oscillation solutions 
predict that the day and the night seasonal dependences will be the same. Therefore, 
in principle one could distinguish between vacuum oscillations and LMA oscillations by 
comparing the seasonal dependence observed during daytime with the seasonal dependence 
observed at night. 

The LMA solution predicts for experiments measuring the low energy 7 Be neutrino line 
(BOREXINO, KamLAND, LENS) that there will not be a significant contribution to the 
seasonal effect beyond that expected from the eccentricity of the earth's orbit. This is 



because the day-night difference is very small at low energies in the LMA solution fl28| , |12 



For vacuum oscillations, there can be significant seasonal dependences of the 7 Be line in 
addition to the purely geometrical effect . 

SuperKamiokande data can be used to test for an enhancement of the seasonal variations 
as the energy threshold for selecting events is increased. Such an enhancement might be 
produced by vacuum oscillations. Vacuum oscillations with Am 2 ~ 4 x 10~ 10 eV 2 , which 
give the best description of the recoil energy spectrum, provide a distinct pattern for the 
enhancement: the effect of eccentricity is larger than the effect of oscillations for a threshold 
energy of 6.5 MeV, while the effect of oscillations is comparable with the eccentricity effect 
for a threshold of 11.5 MeV g§. 
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Comparing Fig. |H] and Fig. || we see that the enhancement with energy threshold 
predicted by LMA is weaker than for vacuum oscillations. For LMA, the increase in the 
threshold only enhances the seasonal variations by of order 30% to 40% of an already small 
effect. This insensitivity to changes in energy threshold is in agreement with our calculations 
of the day and night recoil energy spectra (see Sec. |1V B| ). 

Practically speaking, LMA predicts for SuperKamiokande that there will be no measur- 
able change in the seasonal effect with increasing energy threshold. 



VI. SUMMARY AND DISCUSSION 
A. The current situation 

The Large Mixing Angle solution is consistent with all the available data from solar 
neutrino experiments. Figure |l] shows the allowed region of the LMA parameters in the 
approximation of two oscillating neutrinos. 

We have investigated in this paper the LMA predictions for SuperKamiokande of the 
distortion of the electron recoil energy spectrum, the integrated day-night effect, the zenith- 
angle dependence of the event rate, and the seasonal dependences. We have also evaluated 
and discussed an independent test, which we call the day-night spectrum test. Finally, we 
have analyzed the possibilities for distinguishing between LMA and vacuum oscillations (see 
Table [TTJ and the discussion in Sec. |Vp. We previously showed that the LMA solution is 
globally consistent with the measured rates in the chlorine, Kamiokande, SAGE, GALLEX, 
and SuperKamiokande experiments fl2H . 

The electron recoil energy spectrum predicted by the LMA solution is practically undis- 
torted at energies above 7 MeV, i.e., the v e survival probability is essentially independent 
of energy at high energies. However, the recoil energy spectrum that is measured by Su- 
perKamiokande suggests an increased rate, relative to the standard recoil energy spectrum, 
for energies above 13 MeV |§ggg|. 

We have discussed in Sec. [II B| two experimental possibilities for explaining the upturn 



at large energies of the recoil energy spectrum: (1) a statistical fluctuation that will go away 
as the data base is increased; and (2) an error in the absolute energy calibration. Figure |2| 
shows that a best-fit to the upturn at large energies would require a several sigma error 
in the absolute energy calibration. Both of these seemingly unlikely possibilities, a large 
statistical fluctuation and a large error in the energy calibration, will be tested by future 
measurements with SuperKamiokande. 

There are at least two possibilities for explaining the spectral upturn that do not involve 
experimental errors or statistical fluctuations: (1) a hep flux that is approximately 10 to 
30 times larger than the conventional nuclear physics estimate; and (2) vacuum oscillations. 
We have concentrated in this paper on the explanation that invokes a larger-than-standard 
hep flux. Our results are illustrated in Fig. |3] Even when the enhanced rate is used in 
the calculations, the hep neutrinos are so rare that they do not effect anything measurable 
except the energy spectrum. 

At energies below 7 MeV, the LMA solution predicts a slight upturn in the recoil energy 
spectrum relative to the standard model energy spectrum. This effect is shown in Fig. |3| 
and may be detectable in the future with much improved statistics. 
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However, Fig. shows that even a modest error in the absolute energy calibration could 
give rise to an apparently significant energy distortion. In future analyses, it will be impor- 
tant to include the absolute energy calibration error as one of the parameters that is allowed 
to vary in determining the best fit and the uncertainty in the spectrum shape. 

The zenith-angle dependence and the integrated day-night effect observed by Su- 
per Kamiokande are consistent with the predictions of the LMA and, in fact, show a hint 
of an effect (~ 1.6a) that is predicted by LMA oscillations (see Ref. [I5| and Eq. Fig- 
ures [| and § show that the predicted and the observed nighttime enhancement are both 
relatively flat, approximately independent of the zenith angle. Moreover, all five of the mea- 
sured nighttime rates are above the average daytime rate. However, the results are not very 
significant statistically. If the LMA description is correct, then another 5 to 10 years of Su- 
perKamiokande data taking will be required in order to reveal a several standard deviation 
effect. 

The difference between the daytime and the nighttime recoil energy spectra may be 
detectable in the future. Figure |6] shows separately the predicted day and the predicted 
night spectral energy distributions. 

The most useful way to test for differences between the shapes of the daytime and the 
nighttime energy spectra is to normalize both spectra to the same value. Figure [F] shows the 
predicted difference between the spectra observed during the day and the spectra observed 
at night. It will be extremely interesting to test the hint that a day- night effect has been 
observed with SuperKamiokande by comparing, as in Fig. [7|, the equally-normalized day and 
night spectra. The simplest way of performing this test would be to compare the total area 
of the energy spectrum that is observed below 10 MeV at night with the total area observed 
below 10 MeV during the day. The predicted average difference in the day-night spectrum 
test is about 1.5% for Am 2 = 2 x 10~ 5 eV 2 and decreases to ~ 0.5% for Am 2 = 8 x 10~ 5 eV 2 . 

If the LMA solution is correct, then when the day and night spectra are normalized so 
that they have equal total areas, then the area under the daytime spectrum curve below 10 
MeV must be larger than the area of the nighttime spectrum below 10 MeV. No significant 
difference between the daytime and the nighttime spectra is expected if vacuum oscillations 
are the correct solution of the solar neutrino problems. 

Seasonal differences are predicted to be small for the LMA solution. Equation (|10D shows 
that seasonal differences are expected to be reduced relative to the day-night asymmetry 
by a factor of order of 6 for the location of SuperKamiokande, i.e., the difference between 
the average length of the night in the winter and the summer divided by 24 hours. The 
predicted and the observed seasonal dependences are shown in Fig. § and Fig. |9|. If the 
LMA solution is correct, it will require many years of SuperKamiokande measurements in 
order to detect a statistically significant seasonal dependence. 



B. How will it all end? 

If there is new physics in the neutrino sector, then experimentalists need to provide two 
different levels of evidence in order to "solve" the solar neutrino problems. First, measure- 
ments must be made that are inconsistent with standard electroweak theory at more than 
the 3<T level of significance. Second, diagnostic measurements must be made that uniquely 
select a single non-standard neutrino solution. 
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Where are we in this program? 

We are much of the way toward completing the first requirement in a global sense. A 
number of authors have shown |T^|37| -|39| that an arbitrary linear combination of fluxes 
from different solar nuclear reactions, each with an undistorted neutrino energy spectrum, is 
inconsistent at about 3a or more with a global description of all of the available solar neutrino 
data. The data sets used in these calculations have gradually been expanded to include the 
results of all five solar neutrino experiments: Homestake, Kamiokande, SAGE, GALLEX, 
and SuperKamiokande. The results have become stronger as more data were added. The 
precise agreement |T3| of the calculated sound speeds of the standard solar model with 
the measured helioseismological data has provided a further argument in support of the 
conclusion that some new physics is occurring. 

However, we do not yet have a measurement of a smoking-gun phenomenon, seen in a 
single experiment, that is by itself significant at a many sigma level of significance. This has 
been the goal of the current generation of solar neutrino experiments. 

What can we say about the possibility of achieving this goal with SuperKamiokande? 
There may of course be new physical phenomena that have not been anticipated theoretically 
and which will show up with a high level of significance in the SuperKamiokande experiment. 
We cannot say anything about this possibility. 

If the LMA description of solar neutrinos is correct, we can use the results shown in 
Fig. 0-Fig. HO to estimate how long SuperKamiokande must be operated in order to reduce 
the errors so that a single effect (spectrum distortion, zenith-angle dependence, or seasonal 
dependence) is significant at more than the 3a level. Figures |^[10| show that the errors must 
be reduced for any one phenomenon by at least a factor of two in order to reach a multi- 
sigma level of significance. Since the available data comprise 708 days of operation, over a 
calendar period of order three years, it seems likely that SuperKamiokande will require of 
order a decade of running in order to isolate a single direct proof of solar neutrino oscillations, 
provided the LMA description is correct. 

Fortunately, SuperKamiokande can make a global test of the standard electroweak hy- 
pothesis that nothing happens to neutrinos after they are created in the center of the sun. 
This global test could become significant at the several sigma level within a few years even 
if LMA is correct. For example, the combined effect of the zenith-angle dependence plus 
a possible spectral distortion at low energies might show up as a clear signal. Or, both 
the integrated day-night effect and the day-night spectrum test could be observed. The 
combined statistical significance of several difference tests could be very powerful. With the 
great precision and the high statistical significance of the SuperKamiokande data, we think 
that there is a good chance that a many sigma result may be obtainable in a relatively few 
years. 

Note added in proof. Recently, the SuperKamiokande collaboration has made available 
the data for 825 days of observations (T. Kajita, talk given at the Second Int. Conf. "Beyond 
the Desert," Castle Ringberg, Tegernsee, Germany, June 6 - 12, 1999). In this more complete 
data set, the significance of the excess of events at high energies has further decreased 
and the Day-Night asymmetry has increased to approximately the 2a level. Both of these 
developments strengthen slightly the case for the LMA solution. 
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FIG. 1. The allowed region of the LMA MSW parameter space. When only the rates in 
the chlorine, SuperKamiokande, SAGE, and GALLEX experiments are considered, the vertical 
continuous lines are part of the allowed contours at 90%, 95%, and 99% C.L. The dashed contour 
corresponds to the allowed region at 99% C.L. when both the total rates and the Night-Day 
asymmetry are included. The best fit parameters, indicated by a dot in the figure, are for the 
rates-only fit: sin 2 (2(9) = 0.72 and Am 2 = 1.7 x 10~ 5 eV 2 . The best-fit parameters for the combined 
fit are sin 2 (2#) = 0.76 and Am 2 = 2.7 x 10~ 5 eV 2 . The approximately horizontal lines show the 
contours for different Night-Day asymmetries (numerical values indicated); see Ref. [|30| for an 
illuminating discussion of the Night-Day asymmetry. The data are from Refs. [|l]j3|,f||J. 
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FIG. 2. The effect of a hypothetical energy calibration error on the shape of the electron recoil 
energy spectrum. The data points represent 708 days of SuperKamiokande observations ||. The 
three curves were calculated from Eq. (|3|) for constant values of the energy offset parameter, 5, of 
80 keV, 160 keV, and 320 keV. For visual convenience, each curve was calculated with a somewhat 
different normalization. 
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FIG. 3. Distortion of the 8 B electron recoil electron spectrum for large mixing angle MSW 
oscillations. The figure shows the ratio, R, of the measured number of electrons to the number 
expected on the basis of the standard model |l3| [see Eq. (Q)]. Solutions are shown with the 
standard model fluxes of 8 B and hep neutrinos (the approximately horizontal lines) and with the 
standard hep flux multiplied (see Ref. by factors of 8, 15, and 30 (ratio increases at highest 
energies). The calculated ratios are shown for three different values of Am 2 and with sin 2 26 = 0.79 
in order to illustrate the range of behaviors that result from choosing neutrino parameters. Results 
similar to Fig. |3] are obtained if sin 2 26 is allowed to vary over a representative range of the allowed 
global LMA solutions that are consistent with the average measured event rates. The experimental 
points show the SuperKamiokande results for 708 days of observations, Ref. [fj. 
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FIG. 4. The predicted versus the observed zenith angle dependence of the total event rate 
above 6.5 MeV. Here cos @n = cos(-7r — 0), where ©at is the nadir angle and is the zenith angle. 
The first bin is the average of the daytime rate of all of the observations. The figure shows the 
dependence of the v — e scattering rate in the SuperKamiokande detector upon the nadir angle 
of the sun at the time of the observation is shown for a representative range of Am? and fixed 
sin 2 26. The results are averaged over an entire year. The data are from Ref. ||. The calculational 
procedures are the same as in Ref. [28]. 
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FIG. 5. For a fixed value of Am 2 and a range of values of sin 2 26, the predicted versus the 
observed zenith angle dependence of the total event rate above 6.5 MeV. Same as Fig. |] except for 
different oscillation parameters. This figure makes clear the difficulty of reducing the error bars so 
that they are smaller than the size of the predicted effect. 
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FIG. 6. The predicted day and night electron recoil energy spectra for a representative range 
of Am 2 and sin 2 20. The night spectra are indicated on the figure by N; the corresponding day 
spectra are not marked. For the upper panel, the values of Am 2 are (in units of 10~ 5 eV 2 ) 7.9 
(thick solid line), 4.0 (dashed line), and 1.6 (thin solid line). For the lower panel, the values of 
sin 2 20 are 0.69 (thick solid line), 0.87 (dashed line), and 0.95 (thin solid line). 
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FIG. 7. The Day-Night Spectrum Test. The figure shows the characteristic difference between 
the recoil energy spectrum that is calculated for daytime observations and the energy spectrum 
that is calculated for nighttime observations. The particular curve that is shown in the figure was 
calculated for Am 2 = 2 x 10 _5 eV 2 and sin 2 26 = 0.8; similarly shaped curves are obtained for other 
LMA solutions with the same Am 2 but different values of sin 2 29. The difference between the day 
and the night spectral shapes decreases as Am 2 increases. 
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FIG. 8. The predicted seasonal dependence of the total event rate with eccentricity removed. 
The dependence of the v — e scattering rate in the SuperKamiokande detector upon the season of 
the year is shown for a representative range of LMA solutions with variable Am 2 and sin 2 20 = 0.8. 
In order to isolate the seasonal effect of neutrino oscillations, the effect of the eccentricity of the 
earth's orbit has been removed. Figure S refers to all events above 6.5 MeV. 
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FIG. 9. The predicted seasonal plus eccentricity dependences of the total event rate. This 
figure is the same as Fig. || except that for the present figure we have not removed the effect of the 
eccentricity of the earth's orbit in the predictions or the observations. The data are from Ref. ||. 
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FIG. 10. The predicted seasonal plus eccentricity dependences of the total event rate for 
= 11.5 MeV. This figure is the same as Fig. |9| except that the threshold recoil electron energy 
is larger. 
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